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ABSTRACT 

Based upon the formalism of conformal field theory with a boundary, we give a 
description of the boundary effect on fully developed two dimensional turbulence. 
Exact one and two point velocity correlation functions and energy power spectrum 
confined in the upper half plane are obtained using the image method. This result 
enables us to address the infrared problem of the theory of conformal turbulence. 
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Turbulence arises in many physical systems, which obey simple but nonlinear 
equations of motions. 111 The statistical nature of turbulence together with its non- 
linear, nonequilibrium structures seem to require certain non-unitary quantum field 

[2] 

theories for the theory of turbulence. The appearance of divergences in conven- 

[3] 

sional perturbative approach suggest the use of renormalization group technique. 
However, so far this has been unsuccessful leaving turbulence still an outstanding 
problem in theoretical physics. 

Recently, Polyakov 14 ' ' has proposed a new approach based on the assumption 
of conformal invariance in the inertial range of fully developed two dimensional 
turbulence. He argued that certain non-unitary conformal field theories (CFT's) 
provide exact solutions of the fundamental equations of turbulence, which are 
consistent with the known energy spectrum.' 61 However, Polyakov's solutions in 
terms of correlators of non-unitary CFT are plagued with infrared(IR) divergences 
and also highly dependent on the arbitrary vacuum expectation values of operators 
in the theory. Instead of taking the IR problem as a weakness of his approach, 
Polyakov 14 ' 51 has conjectured that the resolution of the problem is actually the most 
important part of the theory. 

In this letter, we introduce a physical boundary to a turbulent system in or- 
der to regulate the IR divergence. Using an 'image method', we determine exact 
one and two point velocity correlation functions and energy power spectrum which 
depend on the distance from the boundary. Away from the boundary, these func- 
tions quickly become isotropic which is in accordance with the Kolmogorov's local 

[7 2] 

isotropy assumption. ' This allows us to use the boundary not only for handling 
the IR divergences, but also for an effective tool to understand the IR behavior of 
the theory. We attempt to give a physical interpretation to these results. 

In two dimensions, the fluid motion is governed by the Navier-Stokes equation: 
u + tapdaipdpd 2 ^ = ud 2 uj + stirring force, (a, (3 — 1,2) (1) 
where uj and ip denote the vorticity and the stream function respectively and v 
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is the viscosity. Stream function is related to the vorticity and velocity through 
u = d 2 ip and v a = e a sdstp. A typical stirring force is provided by letting fluid 
flow through a grid consisting of very thin bars. For fluid with a sufficiently large 
Reynolds number, fully developed turbulence emerges 30 to 40 times the grid bar 

[2] 

spacing away, where the stirring force is absent. The boundary which we use in 
this paper is parallel to the direction of the incoming flow. 

In the following, we first review the case without boundary. Polyakov's CFT 
approach is to identify the stream function ip with a primary field of a certain 
minimal CFT. In the static and inviscid case, the Navier-Stokes equation reduces 
to iL e a pd a %l)dpd 2 il) n = 0. This equation has meaning only when we exercise care 
in defining the operator product at the coinciding points, e.g., by using the point 
splitting method, 

u e af 3d a ilj(z)d l3 d 2 ilj(z) n = lim e a pd a ^{z + a)dpd 2 ip(z), (2) 

where lim implies angle averaging. In CFT, the operator product expansion (OPE) 
of ip has the following structure: 

ip(z + a)i){z) = (aa) A *~ 2Ai ' {(j)(z) + descendents} + • • • , (3) 

where <fi is the minimal dimension operator of dimension which is most relevant 
as a — > 0. This leads to 

u e ap d a ^(z)d p d 2 ^(zy ~ lim(aa) A ^ 2A ^L_ 2 I 2 1 -I_ 2 L 2 1 ]0(^), (4) 

with L_ n being Virasoro generators. The right hand side of Eq.(4) vanishes, i.e., 
the Navier-Stokes equation is satisfied, for the following two cases: First is when 
[L-2L_i — L-2L 2 _ i\<p= 0, which for example is satified when is degenerate on 
the level two. The simplest concrete example of this is the minimal model of 
p = 2, q = 5, which we denote Ai(2,S), the critical Yang-Lee edge singularity. 
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Second case is when A^ > 2A^. This together with a further restriction coming 
from the constant enstrophy flux condition 191 has led Polyakov and others find series 
of exact solutions of conformal turbulence.' 4 ' 5 ' 10 131 The first case of (2,5) does not 
have constant enstrophy flux, thus requires a subtle balance between vorticity and 
energy injection. 

In order to calculate the energy spectrum E(k), we need the two point function 

of i>\ 

w^>~1^ F 7+ k - < ff-"« + -- (5) 

which results in the power law behavior 

E(k) ~ k a , (6) 

where a = 4A^ + 1 for (0) = 0, and a = 4A^ - 2A^ + 1 for (0) ^ 0. Unlike in 
unitary cases, one point function (<p) is not necessarily zero in non-unitary theories, 
and as we will see in the following, (0) can be determined by boundary condition. 

In this letter, we will focus on -M(2,5) minimal model restricted to the upper 
half plane (y > 0) to investigate boundary effects on fully developed turbulence 
explicitly. One way to see the effect of a boundary on a turbulent flow would be 
to imagine that each eddy in turbulent flow has an image eddy on the other side 
of the boundary, similar to image charges of electrostatics near a conductor. This 
is consistent with the fact that vertical component of fluid velocity vector vanishes 
at the boundary. Furthermore, this picture enables us to extend Polyakov's CFT 
approach to turbulence to the case with boundary. 

Cardy' 141 introduced CFT with boundary utilizing the image method to de- 
scribe surface critical behavior in two dimensions. He has shown that the restricted 
conformal transformation which preserves boundary conditions can be achieved by 
extending the energy momentum tensor T(z) analytically in the lower half plane 
through T(z) = T(z), when y = Im z < 0. One consequence of it is that for the 
degenerate conformal field theories, the n-point function {<f>(z\, zi) • • • <f)(z n , ~z n ))b 
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can be obtained systematically from the bulk 2n-point function, where (•••)& de- 
notes correlation functions in the upper half plane. This is so because the n-point 
function {<f>(zi,zi) ■ ■ ■ <f>(z n , ~z n ))b satisfies the same differential equation as does the 
bulk 2n point function consisting of charges in the upper half plane as well as their 
images in the lower half plane. 

(2,5) minimal model has only one nontrivial primary field 0(1,2)) °f conformal 
dimension A( 12 ) = —1/5, whose normalization is such that (0(i ] 2)(- 2 l)0(l,2)(- 2 2)) = 
1-21 — ^2 1 4//5 . Here we identify % times the primary field 0(1,2) as the stream function 
■0. The factor % is required for the following reason; if we neglect the factor i, we 
can not simultaneously have a real one point function and a positive two point 
function of velocity, which are necessary in order to be physical. 

Following the Cardy's image method, we obtain the one point function of ip 
with boundary from a two point function of i]){z) and its image ip(z) across the 
boundary such that 

(ip(z, z)) b = {ip(x, y)) b = d^y 2/5 , (z = x + iy,y > 0), (7) 

where is an arbitrary real constant. In the unitary case with boundary, 
is fixed by the asymptotic behavior of two point fuction: (ip(xi,yi)ip(x2,y2))b ~^ 
('4 ; (yi))b(' l P(y2))b as \ x i~ x 2\ — > 00, whereas in the non-unitary case d^ is in general 
arbitrary. The mean velocity can be obtained from the one point function of tp 
such that 

2 

(v x (x, y))b = (d y ip(x, y)) b = -d^y~ 3/5 , ^ 
(v y (x,y)) b = -(d x ijj(x,y)) b = 0. 

We interpret this result as describing a system with a net drift velocity profile of 
turbulent flow parallel to the boundary. The drift velocity in Eq. (8) goes to infinity 
as y — > 0, so it should be rounded off over a certain distance from the boundary. 

In order to calculate the energy spectrum of the system, we need the two 
point function of velocity. Again, this can be obtained from the bulk four point 
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function of ip involving two in the upper half plane and two images in the lower 
half. For Ai (2,5)1 which is degenerate on the level two, the bulk four point function 
satisfies a second order differential equation. Solutions of the differential equation 



can be put in terms of linear combinations of hypergeometric functions of variable 



\Zl-Z 2 \ 
%l2/2 



(^(2l,2l)^(22,22))l 



(- 



^13^24 



V 212223^34^14 



-2/5 



(aF 1 (0 + bF 2 (0), 



(9) 



where we use the notation 23 



21, 24 = 22, Zij = Zi — zj, and F\ 
F2 = £ -1 / 5 2^l(f , §; |; 0- a an d & are constants which can 1 
looking at the limiting behavior of the correlator as \z\ — 22 1 — 
the right hand side of Eq. (9) becomes 



cw 3 4.6. 
2-^H5' 5' 5'' 



;0, 

by 

0. In this limit, 



(^(2l, 2l)^(22, 2 2 ))& 



Zl - ^2 



4/5 



-1/5- 



(10) 



where y — (yi + 2/2) /2. On the other hand, the OPE of tp(zi)tp(z2) as 21 — > 22 
gives 



ipipi 



\Z\ - 2 2 



4A, 



a 



\Z\ - 2 2 



4A^-2A. 



- 1 21 - 2 2 1 4/5 + 1 21 - 2 2 | 2/5 C^ci V) l/ 2/5 . 

(11) 

C^tyi and C^^, are structure functions which are determined from the three point 
function of the bulk theory. Explicit calculation shows that C^j = —1 and 
= |C|V2, where C = - r{ %%^{^ /5) Comparing Eq.(10) and Eq.(ll), 
we get a = -1 and b = 2- 2 / 5 d^\C\ 1 / 2 . 

The correlators of velocity components can be obtained by differentiating the 
stream function correlator, for example, (v y {z\)v y {z2))b = d Xl d X2 (ip(zi)ip(z2)) b- In 
order to obtain the energy spectrum E{k) for k 3> ko where ko is the characteristic 
input wavenumber, we consider the form of velocity correlators in the small distance 



6 



limit. When x\ = X2 and y\ = y + e/2, yi = y — e/2 (e<?/), the velocity correlator 
is given by 



(vy(x 1 ,yi)vy{x 2 ,y2))b = V 6/5 (l - Al_ _ + . . 



2 A l^lV2 c -8/5 1( 2/5 A 1 r 7 e4 3 e6 



When the mean velocity is zero (d^ = 0) and when moved far away from the bound- 
ary, i.e., when y >> 1, the isotropic bulk two point velocity correlation function is 
recovered. Similarly, we have 



/ / n / 4 -6/5 / \ 57 e4 243 e6 

Mxi,yiK(x 2 , W )> 6 = -- C + + + . 

25 v " 1 y V 10 y 1800 2000 
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and also (v x (xi,yi)v y (x2,y2))b — 0- It is straighforward to show that the leading 
order terms in Eq.(12) and Eq.(13) are isotropic in the sense that they do not 
depend on the limiting procedure of \z\ — z%\ — > 0. 

The energy density at point (x, y) is given by ^(v^(x,y))i, = J dkE^ xy ^(k). In 
the momentum space, the energy spectrum E^ xy ^{k) is given in terms of velocity 
correlators: 

E M (k) ~ ±J d 2 x'e ik ^-^(v^x',y')v^x,y)) b . (14) 

In the homogeneous case, Et xy \{k) does not depend on the position (x,y) which, 
however in the presence of boundary, is a function of the position. From Eqs.(12)- 
(14) , we get 

E {x , y) (k) ~ vr (Ak 1 !* + Bd^k^) , (15) 

where A = 2- 1 /5i|rg/| &nd B = 2 -3/5^| C |i/2rji/|_ R ig remarkable that 
the boundary effect explicitly shows up through the one point function of the 
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stream function (ip(z))i, = d^y 2 / 5 . In general, when (0) in Eq.(5) does not vanish, 
Polyakov argued that the energy spectrum is described by E{k) ~ /^A^^A^+i 
while for (<f>) b = 0, E(k) ~ k 4A ^ +1 . Above two cases emerge in a unified form 
in our approach, and from this we see that the spectrum near the boundary is 
dominated by E{k) ~ /^A^+i an j ag ^ e distance from the boundary increases, 
the E(k) ~ / c 4A v ,-2A +i behavior becomes more and more dominant. This is 
one of our main result, which seems to be a salient boundary effect in any CFT 
approach to turbulence. It would be very interesting to verify such a phenomena 
experimentally in real turbulent systems. 

One of the basic assumptions in turbulence theory is the Kolmogorov's idea of 
local isotropy, that is, small-scale disturbances in turbulence with sufficiently large 
Reynolds number can be regarded as isotropic, regardless of large scale anisotropy. 
This is consistent with our CFT approach. In our case, large scale anisotropy 
is introduced through the boundary. Nevertheless small scale fluctuations become 
isotropic quickly as we move away from the boundary. In order to see this explicitly, 
we consider the velocity two point function of Eq.(12) rotated by 90 degrees with 
respect to the boundary. This is given by 

+ ^ |c| vva/y/ 5 ( 1 + lJ-|i£; + A£! + ...), 

where the limit is now such that y\ = yi = V and x\ = re 2 + e. Comparing this 
with Eq.(12), we see that the isotropy is broken only at the sixth order of e/y in 
d^p independent terms, and at the second order in d^ dependent terms. 

In this letter, we have determined one and two point velocity correlation func- 
tions by introducing a boundary, which incorporates the large scale structure into 
the scheme of CFT approach to turbulence. In his application of CFT to the tur- 
bulence problem, Polyakov [4,j1 introduced an IR cut-off by adding 5 functions in 
the momentum space. In our approach, a more realistic IR cut-off is introduced 
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(16) 



through a boundary. Although our model defined on a half plane still has an IR 
divergence problem, this can easily be remedied by employing more complicated 
boundaries. Details about other types of boundary and also different minimal 
CFT's leading to more realistic cases will be considered elsewhere. 
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